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, We analyse various phonon-assisted mechanisms which contribute to the decoherence of excitonic 

' qubits in quantum dot systems coupled by the Forster-type transfer process. We show the significant 

04 ' loss of coherence accompanied by dissipation due to charge carrier oscillations between qubit states 

^ , by using a model of one-phonon assisted Forster-type transfer process. We obtain explicit expressions 

^ ■ for the relaxation and dephasing times for excitonic qubits interacting with acoustic phonons via 

' both deformation potential and piezoelectric coupling. We compare the decoherence times of the 

GaAs/AlGaAs material system with half times of concurrence decay in a 2(g)2 bipartite mixed 
excitonic qubit system. We extend calculations to determine the influence of phonon mediated 
interactions on the non-unitary evolution of Berry phase in quantum dot systems. 

PACS numbers: 03.67.Lx, 03.67-a, 78.67.Hc, 73.20.Mf, 85.35.Gv 
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I. INTRODUCTION 



^ - 

The use of excitons (electron-hole correlated states) in quantum dots has generated great interest as suitable devices 
for the solid-state implementation of quantum logic gates^i^i^. The significant advancement in creating and probing 
. excitonic states in quantum dots^i^ and newly developed techniques which allow quantum coupling between excitonic 
states to be continuously varied by external electric fields'' has provided rapid progress in the field of solid-state 
r , quantum computation. In the case of two coupled quantum dot systems, exciton-exciton dipole interactions give rise 
' ^ to diagonal terms which allows quantum logic to be performed via ultra-fast laser pulses on time scales less than 
calculated decoherence times^. Related studies^'^'^ have suggested the importance of the Forster energy transfer 
process in producing entangled excitonic states and using external electric fields to achieve critical logic gate actions. 
Dipole-dipole interaction was first proposed by Forster— and further extended by Dexter— as a mechanism of energy 
transfer in which excitonic wavefunctions localized at different quantum dots do not overlap. 

Other than the actual generation of entangled excitonic qubit states, the successful implementation of quantum logic 
systems depends critically on the decoherence properties of the qubits. Decoherence due to enviromental factors such as 
' phonons and impurities is inevitable and is regarded as a major drawback in solid-state devices. Several strategies such 
[ as decoherence-free subspace a^^d^ ^ optimal control techniqueai^ and immunization processesi^ have been proposed 
' to counter the detrimental effects of decoherence. Geometric quantum computation has gained increased status as a 
robust fault-tolerant scheme in recent yearsi^iii. The geometric phase component (Berry phase) of the phase shift 
in quantum states relies solely on the global geometry of the path executed during cyclical evolution and hence is 
manifestly gauge invariant and consequently insensitive to stochastic operation errorsiS'. It is not immediately clear 
however whether these proposals invoke other routes which facilitate decoherence and if so the role played by various 
system parameters in minimizing errors during logic gate operations. Hence there is a need to scrutinize the various 
intrinsic processes and in particular phonon induced types in order to achieve fault tolerant quantum switches. 
The interaction of charge carriers with acoustic phonons arises mainly from deformation-potential and piezoelectric 
. coupling in quantum dots^^. Theoretical estimates by Fedichkin at a& show that the error rate due to acoustic 
O ' phonons may be a major factor limiting qubit performance, in agreement with experimental results^! of a qubit 
system implemented using a double-dot GaAs/AlGaAs system with a two dimensional electron gas. Electron-phonon 
interaction has been shown be significant in GaAs/InAs quantum dots compared to the bulk system of the same 
material leading to the formation of polarons^'^. In another work, Jacak and coworkers have studied the significance 
^ . , of relaxation of electrons surrounded by a cloud of phonons due to the anharmonic interaction of optical phonons 
' with acoustic phonons in GaAs/InAs quantum dot systems22.. A recent study^"' has considered the possibility of 
localization^^ in which an electron state delocalized over two quantum dots is converted into a mixture of two classical- 
like localized states and consequently loses its coherence to the surrounding reservoir of phonons. 

The interaction of qubits with phonons is generally modelled as a first-order process taking into account only the 
coupling of carrier states localized in quantum dots^- in contact with a phonon bath. The interaction of charge 
carriers with phonons as the qubit flips back and forth between the quantum dots has not been well investigated in 
earlier studies related to phonon mediated decoherence in semiconductor systems^i^i^. The loss of coherence due 
to charge carrier movements during oscillations between the different qubit states is unavoidable, and therefore we 
aim to closely examine this critical mode of decoherence in this work. We model this route of degradation in the 
purity and entanglement of qubit states to proceed via a one-phonon assisted oscillation of qubit states leading to 
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relaxation as well as dephasing of the qubit states. In order to obtain explicit results, we consider the specific case 
of excitonic qubits in Forster coupled quantum dots^. An earlier work^ has examined the coupling of the Forster 
coupled quantum dot system to the external environment consisting of a single radiation mode and the associated 
spontaneous emission and decay processes. We are unaware of other detailed studies of phonon mediated interactions 
involving excitonic qubit states in Forster coupled quantum dots. 

It is not immediately clear whether decoherence phenomenon is dominated by relaxation or pure dephasing process 
and the extend to which the outcome is affected by the nature of qubit-phonon coupling during a cycle of logic gate 
operation. The subtle difference between relaxation and pure dephasing of entangled systems has to be emphasized. 
While relaxation brings about shifts in the population of qubits states with one state favored over another, pure 
dephasing process results in changes in the energy difference between qubit states with the population of qubit states 
remaining intact. While both processes lead to decoherence of the entangled system, we aim to examine the various 
parameters affecting the two inherent sources of dissipation in this work. Our approach undertaken here can be 
extended to other designs of coupled dot qubit systems such as the electrically controlled charge qubit in the double 
quantum dot system or the double impurity semiconductor systema^S. 

It is well known that a quantum system exposed to an environment is best described by the reduced density 
matrix for which approximation techniques'^'^ have been developed to describe the quantum state evolution at times 
determined by the environmental noise. Non-Markovian stochastic Schrodinger equations^, Lindblad^ and Redfield 
master equations'^'^ have been independently used to investigate the dynamics of open quantum system in contact 
with noisy environments. In a recent work, Braur^^ showed that two qubits not interacting directly but which are in 
contact with a common environment can become entangled and the entanglement may last for a long time depending 
the strength of coupling and environmental conditions. These results are similar to those obtained by Ficek et al^ 
who studied a system of two suitably prepared atoms immersed in a vacuum field. The term "sudden death" has been 
coined in an earlier work"^^. to describe the process in which two initially entangled atoms which become noninteracting 
later become completely disentangled in relatively short times. The presence of such features have not been thoroughly 
investigated in excitonic qubit systems. 

In the second part of this work, we investigate the dynamical behaviour of entanglement of excitons coupled to two 
phonon baths associated with deformation-potential and piezoelectric coupling in the absence of a direct channel of 
interaction, i.e the Forster interaction is switched off by control of external electric fields. We follow the approach 
of Tolkunov et al'^^ in which the degree of loss of coherence is determined by the decay of of the off-diagonal terms 
of the reduced density matrix. We determine two important quantities, namely the Berry phase and concurrence^ 
and show the strong correlation between the Berry phase correction fraction which quantify departure from unitary 
evolution and half times of concurrence decay in GaAs/AlGaAs quantum dots. 

II. EXCITONIC QUBITS IN QUANTUM DOTS 

We consider two excitons in their ground states in adjacent coupled quantum dots located at Ra and Rb. We 
assume the quantum dots to be shaped in the form of either cuboid boxes or quasi-two dimensional disks in which 
the vertical confinement energies of charge carriers are larger than their lateral confinement energies. We label the 
localized excitonic states as |Ra) and |Rb). Employing some of the notations in our earlier works^, we represent |Ra) 
as 

1^^^ " ^ 51 *(R-a,re||,rh||,2;e,Z/i) air, a„,rh (1) 
re.rh 

where i/q is the volume of the unit cell, L is the quantization length and aj ^.^ (a^ ) is the creation (annihilation) 
operator of an electron in the conduction (valence) band, denoted by c (u). |0) in Eq. [1] denotes the electronic 
state of the quantum dot in which all electronic ground states are occupied and all excited states are unoccupied. 
For simplicity, we ignore the role of spin effects and spin indices in Eq. [TJ The position vectors Tg and rh are 
decomposed into components parallel and perpendicular to the lateral direction of the quantum dot as = {r^w , Ze) 
and Th = {Yh\\,Zh)- The vertical confinement energies of charge carriers are stronger than the lateral confinement 
energies, the exciton wavefunction and therefore we factorise $(Ra, re|| , rh|| , ^e, ^h) as 

<I>(Ra,re||,rh||,2e,Zh) = ^'(Ra, re|| , rh|| ) (pe{Ze) fh{zh), (2) 

where (/3e(ze) {iph{zh)) is the envelope function of the electron (hole) in the vertical direction of the quantum dot. 
The form of the in-plane exciton wavefunction ^'(Ra, , rh||) at Ra depends on the degree of confinement of the 
electron- hole within the quantum dot. We consider a strong individual charge carrier regime in which the kinetic 
motions of the electron and holes are quantized separately so that the resulting discrete energy levels are affected 
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by the Coulomb interaction between the electron and hole only by a small amount of the order « where L is the 
quantum dot radius. 

4'(Ra,re||,rh||) is further factorized using a suitable change of coordinates 

R - j^(lhreli -Ici-hll) 



where Ir = \/lc^ + Ih^- The effective lengths Ic and Ih are related to the respective electron and hole effective masses 



and confining potential frequencies, and by lo = and Ih = ^^^^ . We obtain ^'(Ra, rejj, rh||) 

S(R, Ra) V'(i"||) where the confining potential in the lateral direction are modelled using harmonic potentials 



1 . 



S(R,Ra) = ^i— exp(--^|R||-Ra|') (3) 
V^^^R 2Lr 



where Lr^ = lo^lh^/lr^- The length scales Ic and Ih associated with the harmonic potential are smaller than the 
effective Bohr radius of the exciton due to the strong confinement of charge carriers in the quantum dot. 

The Gaussian form of function S(R, Ra) in Eq. [3] yields an explicit expression for its two-dimensional Fourier 
transform: 



S(R,Ra) = j(fci\\ exp(iqii-R) S(R,Ra), 



S(R, Ra) - ^ exp(-iq|| . Ra - (4) 

The confining potential in the vertical direction are modelled using harmonic potentials defined by the vertical 
confinement dimensions Izc and Izh respectively for electrons and holes 

1 1 



' V^lzc 21z 



^,(z,) = (^)i/2exp(-^), (5) 

Vl'izh 21zh 

The excitonic state |Rb) is analogous in form to Eq. [T] 

We code the excitonic qubits states using the relative position of the exciton via the basis set, (|L) , |R)) 

|L) = |Ra)®|0), 

|R) = |0)^®|Rb), (6) 

The states, |0)^ and |0)^, which correspond to the absence of excitons denote the respective ground states of the 
quantum dots at |Ra) and |Rb)- We simplify the approach by considering a two- level system involving only the states 
|L) and |R) and work in the limit of a pure Forster coupling. The direct Coulomb interaction which causes the 
formation of the biexciton state |Ra) |Rb) is neglected. We also exclude the possibility of entangled states involving 
the vacuum state, |0)^ |0)b- The two level excitonic qubit Hamiltonian is written as 

Ail 

Hcx~qh ^ -H— CTz + F a^), (7) 

where the Pauli matrices = |L) (R| + |R) (L| and ct^ = |L) (L| — |R) (R|. Aft ~ ila — ^h is the difference in exciton 
creation energy between the quantum dot at Ra and that at Rb- F denotes the interdot Forster interaction amplitude 
responsible for the transfer of an exciton from one quantum dot to the other without involving a tunnelling process. 
The eigenstates of the interacting qubit system in Eq. [7] appear as 

\Xs) = cos(/3/2) |L) + sin(/3/2) |R) (8) 
IXas) = sin(/3/2)|L) - cos(/3/2) |R) , 
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with eigenenergies 



E„ = n„ + Si,-^-^/(^)' + F2, (9) 



r^o denotes the ground state energy of the system in which both the quantum dots are not occupied by excitons 
while Qi denotes a higher energy level of the excitonic qubit system. The energy difference between the eigenstates 



is Y Ail^ + (2F)^ and can be studied as a function of time. In the absence of any decoherence process, the excitonic 
qubit oscillates coherently between the two dots with the Rabi frequency Eg — -Eas- The polar angle f3 in the Bloch 
sphere representation of a qubit is related to Af2 and F via 

2F 

tan/3 =^ (10) 

III. INTERDOT FORSTER AMPLITUDE TERM F 

The magnitude of the interdot Forster amplitude F is given by the matrix element 

F - (R|i?f |L) (11) 
= (5,Ra;/,Rb|C/(|ra-rb|)|/,Ra;5,Rb) 



J2 (.9, Ra; /, Rbl Ui\r, - rtl) Iff, Rb; /, Ra 



where 



C^(lra-rbl) 



e2 



e|ra - rbl 



and Ta and ra are position vectors of the charge carriers with origins at the center of the quantum dots at Ra and 
Rb respectively. The background dielectric constant, e = eq Cr where the relative permittivity, is assumed to be 
independent of the location of charge carriers. In the two-particle interaction matrix element 

(g, Ra; /, Rbl U{\r, - rb|) |/, Ra; 5, Rb) , 

the states to the right of the scattering potential J7(|ra — rb|) represent the initial states while those to the left represent 
the final scattered states. The first matrix term in Eq. [TT]is due to Coulomb interaction in which the excited electron 
in its initial state / in the quantum dot at Ra and a hole in the gth state in the quantum dot at Rb get scattered 
through the potential U to final states in which the electron and hole remain in the same excited states but their 
positions are exchanged. The net effect can be viewed as a tunnelling process in which the exciton gets physically 
transferred from one quantum dot to another. The second matrix element in Eq. [11] represents the scattering of an 
electron in its initial state / in the quantum dot at Ra and a hole in the gth state in the quantum dot at Rb through 
the potential V such that they remain on their respective quantum dot sites but their excited states are exchanged. 
This second process can be viewed as one in which the excited electron in one quantum dot recombines with the hole 
in the ground state and the liberated energy is transferred to excite the electron-hole pair in a neighbouring quantum 
dot. We assume that the quantum dots are separated by a distance, W = |Ra — Rb| such that W |ra — rb| so that 
tunneling effects are negligible. The interdot Forster energy transfer process is then determined by the second matrix 
element in Eq. [TT] 

Substituting expressions for |Ra) and |Rb) (see Ec^. [ij into Eq. [Til and employing the Forster multipole expansion 
of the Coulomb interaction we obtain: 

^^'^^'^^ = 47r6|Rl-RbP ^ ^ ^(Ra,re||,rh||)vI/*(Rb,r;,i,ri,||) (12) 

X (pe(2e) ^h{zh) VVkA) ^r^,r^ ^r'^y^, 

X [/ia ■ Mb - 3(/^a • Rd)(Mb ■ Rd)] 
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where 

_ (Ra - Rb) 

iRa-Rbr 

and 



Ma(b) = J Uc(r) (ra(b) - r) Uy{r) dr 

where Uc(r) and Uv{r) are the respective periodic components of the bulk Bloch functions of the electron and hole. 
We assume that the Bloch functions are identical in both the quantum dot and barrier materials. 
Using Eqs. [3] -[5] in Eq. [T^l we obtain: 

F(W) = (R|Hf|L) (13) 

TT ^2 / 21^1^ \ / 21,elzh ^ 



8 eW3 Vlc'+lhV Vlzc'+W 
where fi = — fJ'b a-nd W — |Ra — Rb|- From Eq. [H]one obtains 

(XsI^fIXs) = -(xasl^FlXa.) =sin/3F(W) (14) 
{x.s\Hf\Xs) - -cos(3F{W) 

In view of the fact that F is key to the entanglement and the formation of entangled states, |xs) and |xas), we expect 
the inevitable loss of energy to lattice vibrations coupled with the excitation transfer to lead to decoherence and 
subsequent degradation of qubit states. We will examine this important process in greater detail in the next section. 



IV. EXCITONIC QUBIT-PHONON INTERACTION 

We consider the exciton state at each quantum dot to be coupled to a continuum of acoustic phonons via both 
deformation potential and piezolelectric coupling. The acoustic phonons are modelled as bulk modes which is a 
valid approximation for quantum dots fabricated in barrier materials with almost similar lattice properties such as 
(In,Ga)As quantum dots with a (Ga,Al)As barrier. The total Hamiltonian of a system of Forster coupled quantum 
dot interacting with phonons is given by 

^qb™ = -f^cx-qb + -H^'^'' + -ff<S;^ qb (15) 

where Hcx-qh is given by Eq. [7|and the matrix elements of Hp is defined by the Forster amplitude term F in Eq. 1131 
H^^ is the phonon reservoir 

q 

where 6^(q) and b\{q) are the respective creation and annihilation operator of a A— mode phonon with wave vector 
q. The A-mode is denoted LA for longitudinal acoustic phonons and TA for transverse acoustic phonons. The 
acoustic phonon energy spectrum is determined by the dispersion relation o^^^a = "^LAlql ^'^^ ^^e longitudinal mode 
and to — ^TAlql foi' the transverse mode, with and u^,^ denoting the corresponding sound velocities. 

is the Hamiltonian for excitonic qubit-phonon interaction via deformation potential coupling and is linear 
in terms of phonon creation and annihilation operators 



^fx -qb = E \/ o ,1 [Mr + Mp a.] (16) 



X 



(6l(-q) + 6A(q)) Kx){n^x\ 
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where |riqA) denotes the occupation number of A— mode phonon with wave vector, q. V is the crystal volume and 
p is the mass density of the material system. and are the respective deformation potential constants for the 
conduction and valence bands. The terms Mr and Mp are written as 

Mr = (Xas;nqA ± 1| {D.e'^ ''^ ^ D^e'^ ''^) \xs;n^x) , (17) 

and 

Mp = (xs;nqA ± l\{D,e'^''^ -D,e''i''-)\xs;n^x) (18) 
- (xas;nqA ± l|(i?ce''i'"--i?.e'i'->')|xas;nqA) 

The Mr term describes a decoherence process in which the qubit state |xs) relaxes to the qubit state |Xas) due to 
mediation by phonons. The Mp term describes a second type of decoherence process in which a shift occurs in the 
energy difference between the two qubit states resulting in pure dephasing of the entangled system. 

A similar expression as in Eq. [16] can be obtained for the Hamiltonian Hcx'^qh a describing exciton-phonon interac- 
tion via piezoelectric coupling 



-ffex-qb.A = J2 TT^hx o:r^7~ i^x.xqyqz + iy,\qxqz + S,,,xqo,qy) (19) 



A,q 



[Nr CTx + Np a^] (fel(-q) +feA(q)) l^qA) (^qAl 



where the relative permittivity, is assumed to be unaffected by the contribution from strain fields associated with 
acoustic phonon modes, en is the piezoelectric constant and ^. ^ is the unit vector of polarization of the A-phonon 
along the i-direction. Excitonic interactions with phonons due to piezoelectric coupling are highly anisotropic in 
nature^ and the form of H^^™^ x depends on the choice of polarization components and the modes associated with 
A. Further details of the dependency on LA, TAl and TA2 modes are given in the next section. Nr and Np are 
analogous to Mr and Mp given in Eqs. [T71 and [T51 respectivelv 

Nr - (Xa.; nq A ± 1| [e'^ '"'' - e'^-^-) |x.; a) , (20) 

and 

Np - (xs;nqA±l|(e''i-'-^-e''i-'')|x.;riqA) (21) 
- (Xas;nqA±l|(e"'-"^^-e"i->^)|Xas;r.qA) 

Unlike the Mr and Mr terms, both Nr and Np vanish as q ^ due to exact cancellation of electron and hole form 
factors as the piezoelectric coupling constant remains the same for electrons and holes. 

We consider the Hamiltonian in Eq. [12] to be associated with a interdot Forster energy transfer process 
that proceeds via a one-phonon assisted transfer of excitonic qubits. A one-phonon-assisted decoherence of qubit 
states is more likely at lower temperatures while multiphonon phonon processes are expected to dominate at higher 
temperatures. We express as arising from second-order process involving two Hamiltonians 

Ht = ^''x-qb © ^^F (22) 

where X = DP(Piez) for coupling via deformation potential (piezoelectric fields). The relaxation of qubit states occurs 
via several routes 

(23) 



(24) 



Xs; 


riq 


a) 




|Xas;"q A ± 1) 


^ |Xas;"-q A ± 1 


Xs; 


riq 


a) 




^""^ |Xs;»^q A ± 1) ^ 


^ |Xas;?^q A ± 1) 


Xs; 


riq 


a) 


Hf^ 


1 \ cx — qb 

|Xs;'^qA) — > 


|Xas;"-q A ± 1) , 


Xs; 


nq 


A> 




|Xas;nqA) > 


|Xas;»^q A ± 1) 
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Likewise, pure dephasing of excitonic qubit state \xs'i n-q) can occur via 





cx - qb 1 1 1 \ -^F 1 1 1 \ 

' |Xas;nqA±l) > |Xs;"-q A ± 1) , 






IXs; riq A ± 1) ^ Ixs; riq A ± 1) , 




Xs;nqx) 


— ^|Xs;n-qA> |Xs;n-q A ± 1) , 




Xs;nqx) 


^|Xas;n-qA) IXsi^qA ± 1) 


(25) 



Pure dephasing of the qubit state |xas; "-q) occurs similarly as shown above for the |xs; "-q) state. 



V. EVALUATION OF MATRIX ELEMENTS 

Using Eqs. [1] [31 HI and O in Eq. 1161 we obtain explicit forms for the following matrix elements 



(L;nq A ± l|-ffeTqb |L;nqA) = T.D{q\\,qz) 



,-iqRa 



L;nqA±l|ffex-qb|R;nqA> - ^D{q\\,q.) e-'^ —e ^ 



(26) 
(27) 



where Lr is defined below Eq. [3]and a similar expression as in Eq. [Oilcan be obtained for (R; nq a ± 1| -f^ox-qb l-^^! "^qA) 
The function J^Dlq^^^ qz) which is given by 



y 2pVujqX 



112 2 ll 2 2 1, 2„2 ll 2 2 



differs from the commonly used form^*^ as we consider the non-isotropic propagation of phonon here. From Eqs. [8l 
[26l and [27l we obtain 



1 



(Xas;nqA±l|i/°P b|Xs;riqA) - sin(/3/2) (c 



-iq Ra _ -iq-Rb 



^D{q\\,qz) 



(28) 



-qb l/v-si '-11^/ 2 ' 

- cos(/3/2)e ^I]c((7||,9,) 

The first term in Eq. [28] contains the coherence factor (^e~"''^=' — e^'i'^t) which remains effective over large distances 

between the quantum dots. The second term involves the factor e "^i^ which is short-ranged as it depends on the 
overlap integral of the exciton wavefunctions in quantum dots located at different sites. By integrating the coherence 
factor over the polar angle 0, we obtain 



(Xas;"qA ± l|-ffcx-qb IXs;?^qA) I = '^T^^ oiih qz) 



1 _ w2 

-sin2(/3)(l - Jo(g||W)) + cos2(/3)e ^ 



(29) 



where Jo(g||W) is the zeroth-order Bessel function. In order to obtain a tractable form of Eq. [511 we approximate 
Jo{qs\n{9)W) by taking the angular average 



Jo(gsin(0)W) 



/ / Jo(qsin(6')W) s\n{e)ded4> 



J J sin{9)ded4> 



sin(gW) 



For large separations between the quantum dots, we approximate 



(Xas;?^qA ± l|i?„^_ql, |Xs;n-qA) I ~ TtS^, (g|| , ) slu (/3) (1 - 



sin(gW) , 



Similar expressions as in Eq. [511 but with slight variations in the dependency on polar angle [3 can be obtained for 
I (Xas; "-q A ± 1| ^c^-qb IXas; "q a) ? and | {xs] riq A ± 1 1 qb IXs; "-qA) ?■ Thcsc exprcssious win be used in Eqs. [13 
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and [18] to evaluate the rates at which the excitonic qubits lose their decoherence via relaxation and dephasing in the 
following section. 

In the case of qubit-phonon interaction via piezoelectric coupling, we use a choice of phonon polarization components 
with respect to the cubic crystallographic axes of the zinc-blende type crystal. Substituting Eqs. [1] [31 2] and [5| into 
Eq. I19[ we get explicit forms for the following matrix elements 



(L;nqA±l|ffr.%,A|L;nqA> - Sp(g|| , g,, A) e-'ir^^ (30) 
where the function I]p(g|| , g^, A) is defined as 



(R, -|-R,. ) W 

L;nq;,±l|i?P'-,^JR;nqA) = Sp(g|| , g., A) e^'^i-^^e ^ (31) 



,^,A) = A,(0,0) e~i^-V e-i'c^^ e-^'-^^S] (32) 

eoe^ y 2pVvxq L J 

The anisotropy factor k\{9, (j)) is given by 

3 

Ala(6',(/>) = - sin 26* sin 6* sin 20 , (33) 

Atai(6',</>) = ^ (sin 61 -sin 361) sin 20 , (34) 
8 

P^TA2{0,4>) = i sin 26* cos 20 (35) 

The TAl and TA2 modes correspond to the two possible polarized direction of the transverse phonon. Calculations 
of the square amplitude term | (xas; a ± 1| -ff^lqb.A IXs; "qA) P where A = LA, TAl, TA2 are slightly involved due 
to the presence of the polar angle in Eq. [33|-[35l For A ~ LA, we obtain 



(Xas;"-qLA±l| -ffcx-qb,LA IXs;"qLA) P = ^Sp (g|| , ^2 , LA) 
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sin^(/3)(l-^) + 9cos^(/3)e- 



(36) 



where I]p(q|| , g^, LA) can be obtain using Eqs. 1321 and [551 Similar expressions as in Eq. [SSI can be obtained for the 
TAl and TA2 modes of the transverse phonon using Eqs. [Mland[55l 

VI. DECOHERENCE RATES OF EXCITONIC QUBITS 

The decoherence rates associated with relaxation and pure dephasing processes are calculated using 

1 In 2 

— = Y 51 K/l^intN)| (Afq,A + l)<5(A^±fil^qA) (37) 

where \i) and |/) denote the initial and final states which are functions of |xas) and |xs) (see Eqs. [171 [HI [20l and 
[^T]) . AE is the energy difference between the initial and final states which determines the wavevector of the emitted 
phonon. _ffint is the interaction operator which is substituted by -H^^^qb or ff^'lqb a depending on the type of phonon 

coupling. In the case of decoherence associated with Forster transfer process, H\nt is given by as will be detailed 
in the next section. iVq>, is the thermalised average number of phonons at the temperature T, and is given by the 
Bose-Einstein distribution N^^x — exp(/ia;q A/fcsT) where ks is the Boltzmann constant. 

From Eqs. \Vf \ [TBI [29l and [37l we obtain explicit expressions for the rate of relaxation {^/r'^-p) and rate of dephasing 
(1/r^p) of excitonic qubits 

where the functions Sr{P) and Sp{P) are given by 

SAP) ^ I sin^(/3) ( 1 - + cos^(/3)e-* (39) 



2 V gw y 

SAP) 4cos^(^)(l-^)+isin^(/3)e- 
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To simplify the derivation of Eq. [38l we have used the assumptions that lo « 1 



the imaginary error function, Al = q^J^{^P~-\^ and q ~ 
analytical expression for F(W)). It can be noted that the cut-off frequency in Eq. 
approximately the inverse phonon flight time through the quantum dot. 



1 and Izo ~ Izh = Iz- Erfi(A/) is 
7^ where 7 = ^^^^1 (see Eq. [13] for an 



2F(W) 

I occurs at uji 



^/\ which is 
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FIG. 1: Relaxation time r^^ (dashed) and dephasing time, t^^ (full) as functions of 7 = 2f(w) 
Ic ~ Ih = 1 = 2.5 nm and Izc ~ Izh = Iz = 2.5 nm. 



at T = 10 K, W = 5 nm , 




Dot size (nm) 



FIG. 2: Relaxation time, t^^ (dashed) and dephasing time t^^ (full) as functions of dot size 1 at 7 = 0.25, W = 4 nm, Iz = 1 
nm and T = 10 K. 

Using the parameters for the GaAs/AlGaAs material systemSS^, we calculate the relaxation time t^^ and dephasing 
time as functions of 7 in Fig.[TJ At 7 = 0, the qubit states are maximally entangled and decoherence is dominated 
by the relaxation process as one can expect from two equally populated states. As 7 increases, the difference in the 
population of the qubit states is enhanced and decoherence becomes increasingly dominated by pure dephasing process. 
We note that the minimum in relaxation time is reached at 7^ and is given by t'^ « W/2ttv^j^ thus satisfying the 
condition where the interdot separation matches the phonon wavevector. At small values of 7, the quantum dot size 
parameter 1 has a stronger influence on the dephasing process as shown in Fig. [2] where the dephasing time decreases 
sharply after 1 « 2 nm while the relaxation time increases gradually mainly determined by the exponential term in 
Eq. Fig. [3] shows that at critical values of W (depending on 7 and 1), the relaxation times can reach very small 
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FIG. 3: Relaxation time t^^ as function of interdot distance W a,t ^ — 1.5 (dashed) and 7 — 0.1 (full) at 1 = 2 nm, Iz = 1 nm 
and T = 10 K. 



values in the order of several picoseconds. 

In the case of qubit-phonon interaction via piezoelectric coupling, we obtain the relaxation and dephasing rates for 
A = LA, using Eqs. [TTl [H EHl and [37] 



LA 



l-exp(-i(r2-l)g2l2) 



G,.(^^^)(A^q,A + l)^;(rt(/3) (40) 



where r — Ih/lc and we have used the assumption that 1 = lo = hi order to obtain an analytical expression. The 
function (x) is given by 

^ , , 2x + 15 v^e-^(4a;(x + 3) + 15)Erfi(Vi) 

G.Ax) ^2 (41) 

and S'j.{(3) and S'p{(3) are given by 



SUP) ^sin^(/3)(l-5^^) + 9cos^(/3)e-^ (42) 

S;iP) =cos^(/3)(l-^) + ^sin^(/^)e"^ 

Similar expressions as in Eq. |40] but with terms G.j,^j^ (x) and G^^^ {x) associated with TAl and TA2 modes are 
obtained for qubits interacting with transverse acoustic phonons 

3 V¥e'-(2a; + 3)Erfi(Vi) 
G...(x) =^ (43) 



^TA2 (•^) 



y/^e-''{x{x{2x + 5) + 12) + 15)Erfi (^i) _ x{x + 2) + 15 
6x'^/2 3x^ 



As expected 1 / Tp.'-P^ = when Ih — Ic due to piezoelectric coupling being a polar mechanism. 

Fig Hlshows that dephasing process due to qubits interacting with phonons via deformation potential is the dominant 
mechanism of decoherence for larger values of 7 (> 2.5) while the relaxation process dominates at small values of 7. 
As expected, these features are similar to those obtained in the case of phonon coupling via deformation potential. 
These results as well as those obtained for the deformation potential case highlight the critical role of 7 and quantum 
dot parameters in influencing the decoherence properties of excitonic qubit systems. Therefore the quality factor 
which determines the number of charge oscillations that can be resolved within the decoherence time can be selected 
using the quantum dot system configuration based on parameters 1, Iz, r, W and T. It appears that phonon assisted 
decoherence can be suppressed by careful choice of system parameters leading to higher fidelity of logic gate operation. 
In the next section we consider another source of decoherence involving one-phonon assisted Forster transfer process. 
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FIG. 4: Relaxation time, r^.^^ (dashed) and dephasing time r^.^^ (full) as functions of 7 at W = 5 nm, Ic = Izc = 2 nm, r — Ih/lc 
= 5 and r = 10 K. 
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FIG. 5: Relaxation time Tp.^^ as function of interdot distance W at 7 = 1.5 (dashed) and 7 = 0.1 (full) at Ic = Izc = 2 nm, r 
= 5 and T = 10 K. 



VII. DECOHERENCE ASSOCIATED WITH ONE-PHONON ASSISTED FORSTER TRANSFER 

PROCESS 

The decoherence rate, — associated with phonon assisted Forster relaxation process is obtained using Eq. 1371 

X,A 

where Hmt is evaluated via 

\{f\Hf\^)f^^^\J2T:J' (44) 

i=l 
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where X = DP or Piez and ^ are the transition amphtudes associated with qubit relaxation. T^^ and T^^ for 
example are obtained using Eqs. and [Ml 



J.1 


1 


rpi 


1 



-qb 



In the case of dephasing mechanism, we use 



" 8 



(45) 
(46) 

(47) 



The transition amplitudes associated with 



where the transition amplitudes T*^ are evaluated using Eqs. 1251 - 
transitions between the Xas is obtained analogous to Eq. |44l 

By neglecting the interference between terms in Eqs. I44[ wc obtain explicit terms for the rates of relaxation and 
dephasing 



1 / r[p) _ (-De — DvYq^ 
' " 32V^ft(l+72)pi;2 



, 3 2,2NErfi(Al),^^ 
exp(--(z2l2) V _ , 1^cF 



Al 



-(iVq,A + l)^.^(p)(/3) 



where 7 = 2_f(W) ^-'^"^ ^^"^ functions S^{(3) and Sp{(3) are given by 



sin'^iP) + -cos2(/3)(3- 



-cos2/3) + sin(2/3)e' 



■sin2(2/3) -sin''(/3) 



sin(qW) 



(48) 
(49) 

(50) 



isin3(/3) 



sin((7W/2) 
qW ' 



The combined rates of relaxation and dephasing for qubits coupled to the LA, TAl and TA2 modes of the acoustic 
phonons associated with the piezoelectric fields are obtained as 



r(p) 

Picz.F 



16e2(l +72)npt;' 



exp(-g2i^/2) 



1 — cxp( 



(r2 - l)q^X 



2l2 



G..(:Tl¥^)(^q,A + l)<.)(/3) (51) 



^2(r2 + 1)' 



where the function G^{x) is obtained using Eq. HTland associated functions, G^^-^{x) and G^j^^{x) (see Eq. 
the functions 5*^ {(3) and Sp{(3) are given by 



<(/3) 



3i(/3)+52(/3)e-^ 
sin(gW/2) _ 



.93 (/3) 



sin((7W) 



and 



(52) 



where ^^(a:), « = 1, . . . , 4 can be obtained using Eqs. [331 EH and[35l 

Fig. [51 shows that pure dephasing dominates at low values of 7 in a one-phonon assisted Forster transfer process 
while the relaxation process is the preferred mode of decoherence at higher values of 7. These features are in reverse 
of those obtained in Fig. [T] where the phonon mediated process takes place in the vicinity of individual quantum 
dots. At very small values of 7, the decoherence times reach values of the order of several picoseconds showing the 
significant loss of coherence which occurs during charge carrier oscillations between qubit states. Fig. [71 shows that 
relaxation via phonon-assisted Forster transfer process is a more likely route for decoherence than the direct mode 
of relaxation at large interdot separations and small 7 values. The notable oscillations in relaxation times is possibly 
due to matching of the interdot separation with multiples of the phonon wavevector as 7 is increased. In Fig. [51 we 
note the reverse effects on the two types of decoherence mechanisms due to an increase in size parameter 1 

Fig. [HI shows that the one-phonon assisted Forster transfer process associated with piezoelectric fields proceeds 
in the same manner as in the model involving deformation potential. Fig. \TQ\ highlights the range of W for which 
decoherence times can match typical logic gate times of around 10 - 50 ps. In the presence of an external electric 
field, energy for both electrons and holes is reduced with larger dots experiencing bigger shifts. We expect changes 
to the order of at most a magnitude to the results obtained here. 
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FIG. 6: Comparison of relaxation time r^p ^ (dashed) and dephasing time r^p p (full) due to deformation potential coupling 
as function of 7 at = 4 nm (full), 1 = 2 nm, Iz = 1 nm and T = 10 K. 
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FIG. 7: Ratio of relaxation times jP^ as function of 7 at = 4 nm (full), 5 nm (dashed), 6.5 nm (dotted) and 1 = 2 nm, Iz 
= 1 nm and T = 10 K. 



VIII. BERRY PHASE AND CONCURRENCE OF THE COMPOSITE EXCITON SYSTEM 

We evaluate two important quantities, namely the Berry phase and concurrence of the composite exciton system in 
quantum dots. We consider a simplified model in which the Forster interaction is switched off after t= and assume 
that there is no correlation between the phonon bath interacting with the correlated electron-hole pairs at different 
quantum dots. The Hamiltonian in Eq. [15] can be solved"^^ to produce an explicit expression for the reduced density 
matrix of the two qubits system at a later time t. The overall initial density matrix (pt(0)) of the system of qubits 
and phonon reservoir, Pph{0) is obtained as 

pr{0) = p,40 f ® pe.(0)« ® p^;,(0) ® pf^iO) . (53) 

where Pex{0)^ and Pex{0)^ are the initial density matrix states of the two-state qubit system (see Eq. [5]). Ppi^{0) and 
Pphi^) density matrices associated with the phonon reservoir in the quantum dots at |Ra) and |Rb), respectively. 
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FIG. 8: Relaxation time, t^^ ^ (dashed) and dephasing time r^p ^ (full) as functions of dot size 1 at 7 = 0.25, W — 4 nm, 1^ 
= 1 nm and r = 10 K. 
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FIG. 9: Comparison of relaxation time r^.^^ ^ (dashed) and dephasing time Tp.^^ ^ (full) due to piezoelectric field as function 
of 7 at W = 5 nm (full), 1 = 2 nm, r = 5 nm and T = 10 K. 



At t = 0, each phonon reservoir is assumed to be thermalized at temperature T 

p^.(o)=n(i- 



(54) 



where I — L,R. 

We consider an initial state given by |xs) — cos(/3(0)/2) |L) + sin(/3(0)/2) |R) (see Eq. |S]) which depends on the 
interdot Forster interaction amplitude F a.t t — through the angle /3(0). We will show later that the evolution of 
f3 with time is dependent on the dynamics of interaction between the qubit states and phonon bath. In the case of 
identical qubits, the state prit) — \Xs) (Xs| evolves according to 



Prit) 



/O 






COS' 





^(3/2) 



sin(/3(0)) e 




-iAOt 



zit) 





isin(/3(0))e'^"*Z(t) 
sin2(^) 







0/ 



(55) 
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FIG. 10: Relaxation time Tp.^^ ^ as function of interdot distance W at 'y — 1.5 (dashed) and 7 = 0.1 (full) at 
nm and T = 10 K. 



— 2 nm, r = 5 



where Afi is the difference in exciton creation energy between the two quantum dots and Z{t) = e~^['^^'^*^+'^^*^*^l 
provides a measure of the decay of the off-diagonal matrix elements when each qubit is isolated from the other, 
but exposed to its own reservoir of phonon bath. The function G{t) = exp(— 2 J -^J{uj) coth.{huj/2kBT) sm^{ujt/2)) 
where J(w) is the spectral density function which yields information about the interaction of the quantum dot with 
the phonons^- 



Jd{^) = ^d{q\\ ; <lz)'^S{uj — cijq) Deformation Potential 
q 

Jp{u!,X) = Sp(g|| , (/z, A)^(5(a; — Wq) Piezolectric Interaction 



(56) 
(57) 



where :z]j{q^^,qz) and S]p((7j| , (72, LA) are given in Eqs. 1281 and 1321 respectively At small cu, J{lo) 
exponent k determines between the cases of ohmic (fc = 1), sub- (A: < 1 
explicit expressions for Jd{'^) and Jp{uj,1jA) are determined using Eqs. 
k = i and k ~ b respectively for Jd(^) and Jp{ijJ, LA). 
The Berry phase of a mixed state is defined as 



J where the 

and superohmic (fc > 1) couplings'*^. The 
51 and [5^ to be in superohmic forms where 



rj = arg 



E 



^E,{0)E,iT){x, mix, (T))e- So ^*(x. (*) Ix. (t)) 



(58) 



where the system frequency Af2 determines the quasi-cyclic path of the geometric phase with time t varying from 
Ej{T) and Xj{T) refer to the eigenvalues and eigenvectors, respectively of the density matrix in Eq. [55l 



to T = 



AO ■ 



At t = 0, G{t) = 0, Z{t) = 1 and \x+{t)) ^ IXs(O)). Substituting Eq. [59]into Eq. Elwe evaluate t] ^ tiq + t]c where 
rjQ = 7r(l — cos(/3)) is the well known result for a unitary evolution at G{t) = 0. The unitary phase correction term 
77c which appears when G{t) 7^ is evaluated using series expansion of G{t). 

Fig. [m shows the increase of Berry phase correction fraction ^ with temperature in GaAs/AlGaAs quantum dots 
coupled to phonon baths via deformation potential. The geometric phase of the qubit state undergoes a greater 
degree of nonunitary evolution within a superohmic environment at larger temperatures. At T = = K, zero point 



where tan(^^) 



E±{t) 

\x+it)) 

'/3(0) 



i ± i (cos^(/3) 

Cos(M)|R) 4 



f exp(-2G(t))sin2(/3))' 
sin( — ) |L) 



(59) 



exp(— G(t)) cot(^^i^). and the evolution of the phase is determined only by \x+(t)) as £'-(0) = 0. 
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FIG. 11: Berry phase correction fraction ^ as function of polar angle /3 due to decoherence associated with deformation 
potential at W = 5 nm, 1 = 1 nm and T = 120 K (fuU), 50 K (dashed) and T = K (dotted) 
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FIG. 12: Berry phase correction fraction — as function of the dot size due to decoherence associated with deformation potential 
(dashed) and piezoelectric coupling (full) at = 5 nm, /3 = ^ and T — 100 K. 



fluctuations associated with a noisy environment contributes to the non-zero correction in Berry phase. At (3 — mr 
where n is an integer, T = vanishes so that the system stagnates and ^ = at the nodal points. At /3 = n-|, T 
becomes very large and that the system never completes its quasi-cyclic path and there is negligible phase correction. 
At other values of the departure from unitary behavior is clearly dependent on the spectral density function G(t) 
as the Bloch sphere is altered according to G{t). There is decreased effect of decoherence effects on r] as the phase 
angle sweeps the second quarter to complete the first half of cyclic evolution. 

Fig-HHshows the gradual decrease of Berry phase correction and hence progress towards a more unitary evolution at 
larger quantum dot sizes for both types of phonon couplings. The difference between the two curves can be attributed 
to microscopic properties of the spectral density functions J_D(a;) and Jp{uj). In both cases, the density function 
rises rapidly with uj and reach peak values at an optimum ujm before starting to decrease at even higher values of uj. 
Fig. [T3] shows the increase of Berry phase correction with interdot distance W which gives a measure of the degree 
of entanglement at t = 0. It is to be noted that W provides a direct link between concurrence (details below) and 
correction to the Berry phase. At large W, the concurrence of the qubit states at t = is weak and the Berry phase 
is less immune to decoherence effects and undergoes a higher degree of non- unitary evolution at t > 0. 
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The evolution of the entanglement between the subsystems Pex{^)^ and Pex{^)^ is determined by a quantity termed 
concurrence, C{t) which is related to the entanglement of formation^Si. For a pure or mixed state, pt-, of two qubits, 
the spin-flipped state is defined as 

Pt = {(Ty <Xi ay) p^ {ay (g) ay) , (60) 
where ay belongs to the set of Pauli matrices. The concurrence is given by^ 

4 

C(ptW) = max{0, 2maxAi-^Aj} . (61) 

where A^s are the square roots of the eigenvalues of the non-Hcrmitian matrix ptPt with each value of Xi being a 
non-negative real number. 

The eigenvalues of the matrix product pxpr for the density matrix Eq. [55] are 

Ai.2 = isin(/3)(Z(t)±l), (62) 

and A3 4 = 0. Hence the concurrence is obtained as 

C(t) = isin(/3)e-4[G^W+«-Wl (63) 

Figs, [m and [T5l show differences in half times of concurrence decay in GaAs/AlGaAs quantum dots which are 
isolated from each other but each exposed to their own phonon baths via deformation potential and piezoelectric 
coupling. The half time during which half of the initial entanglement at t = is lost is less for qubits interacting 
with phonons via deformation potential compared the qubit-phonon system incorporating piezoelectric coupling. 
These differences can be attributed to functions Goit) and Gp{t). The function Goit) associated with the spectral 
density functions J_d(w) (deformation potential) rises with time and reaches a constant value whereas Gp(t) retains 
its monotonic increase with time. Both figures show an increase of half time with quantum dot size 1 which is due 
to the shorter time needed for the phonons to travel the length of the quantum dot and decrease of spectral density 
function J(w) as the dot size is increased. It is to be noted that in our model wc have considered lattice vibrations at 
the two qubit locations to be uncorrelated so that the bath modes are independent unlike the model incorporating a 
common environment used by Braun^. 



IX. CONCLUSIONS 



We have studied several processes that contribute to the decoherence of excitonic qubits in quantum dot systems 
coupled by the Forster-type transfer process. We show the significant loss of coherence that occurs through charge 
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FIG. 14: Concurrence decay in qubits coupled to uncorrelated phonon baths via deformation potential coupling at 1 = 2 nm 
(dashed), 1 nm (full), 1.5 nm (dotted) and T = 10 K, with the prefactor i sin(/3) suppressed. 
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FIG. 15: Concurrence decay in qubits coupled to uncorrelated phonon baths via piezoelectric coupling at 1 = 3 nm (dashed), 
0.5 nm (full) and T = 10 K, with the prefactor i sin(/?) suppressed. 



carrier oscillations between the different qubit states by using a model of one-phonon assisted Forster-type transfer 
process. It is shown that increasing a tuning factor 7 has opposite effects on the relaxation and dephasing times 
for excitonic qubits interacting with acoustic phonons via both deformation potential and piezoelectric coupling. 
Our results show that entanglement of qubits which are isolated from each other last for comparatively long times 
when interacting with uncorrelated phonons baths via deformation potential compared to the case of piezoelectric 
coupling which become completely disentangled within 20 ps. The results obtained here show that phonon mediated 
interactions exert a sizable influence on the non-unitary evolution of Berry phase in quantum dot systems. Lastly, our 
results emphasise the strong correlation between the Berry phase correction fraction which quantify departure from 
unitary evolution and half times of concurrence decay in GaAs/AlGaAs quantum dots. 
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